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Abstract 

Solutions of boundary value problems for a diffusion equation of fractional 
and variable order in differential and difference settings are studied. It is 
shown that the method of the energy inequalities is applicable to obtaining a 
priori estimates for these problems exactly as in the classical case. The cred- 
ibility of the obtained results is verified by performing numerical calculations 
for a test problem. 
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1. Introduction 

Fractional calculus is used for the description of a large class of physical 
and chemical processes that occur in media with fractal geometry as well as 
in the mathematical modeling of economic and social-biological phenomena 
[11, Q, [s], 3, Isf . In general, a medium in which a process proceeds is not ho- 
mogenous, moreover, its properties may vary in time. Mathematical models 
containing equations with variable order derivatives provide a more accurate 
and realistic description of processes proceeding in such complex media (see 

e.g. BBS)- 

Therefore, the development of numerical and analytical methods of the 
theory of fractional order differential equations is an actual and important 
problem. 
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Numerical methods for solving variable order fractional differential equa- 
tions with various kinds of the variable order fractional derivative have been 
proposed 0, Eol, [111, [H [l3|, Q . 

The positivity of the fractional derivative operator has been proved in [l[ 
and this result allows to obtain a priori estimates for solutions of a large class 
of boundary value problems for the equations containing fractional deriva- 
tives. The authors of the paper |15| have obtained a priori estimate for 
the solution of the Dirichlet boundary value problem of a fractional order 
diffusion equation in terms of a fractional Riemann-Liouville integral. The 
fractional diffusion equation with the regularized fractional derivative has 



been studied, for example, in [16|]. In the papers [17|, ll8|, ll9|, the diffusion 



wave equation with Caputo and Riemann-Liouville fractional derivatives has 
been studied. The difference schemes for boundary value problems for the 
fractional diffusion equation both in one and multidimensional cases have 



been studied 
tained in 




A priori estimates for the difference problems ob- 
by using the maximum principle imply the stability and 
convergence of the considered difference schemes. 

Using the energy inequality method, a priori estimates for the solution 
of the Dirichlet and Robin boundary value problems for the diffusion-wave 
equation with Caputo fractional derivative have been obtained |22]. More 



references on fractional order differential equations, including the diffusion- 



wave equation, can be found, for example, in [23 



2. Boundary value problems in differential setting 

2.1. The Dirichlet boundary value problem 

In rectangle Qt = {{x,t) : < x < Z, < t < T} let us study the 
boundary value problem 



^o/^^^M = ( k{x,t)^ ] -q{x,t)u + f{x,t), < X < 1,0 <t <T, (1) 



d 



dx 



du 



dx 



u{0,t) = 0, u{l,t)=0, 0<t<T, 



u{x, 0) = uo{x), < X < I. 



(2) 
(3) 
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Where < Ci < k{x,t) < C2, q{x,t) > 0, d'^l;^\{x,t) = J^Ur{x,T){t — 
r)~"*^^^(ir/r(l — a{x)) is a Caputo fractional derivative of order a{x) , < 
a{x) < 1 , for all x E (0, 1), a{x) E C(0, T) [Isl, Q- 

Suppose further the existence of a solution u{x,t) E C^'^{Qt) for the 
problem ([I])-®, where (7™'" is the class of functions, continuous together 
with their partial derivatives of the order m with respect to x and order n 
with respect to t on Qt- 

The existence of the solution for the initial boundary value problem of a 



number of fractional order differential equation has been proved [25|, |26|, |27 
Let us prove the following: 

Lemma 1. For any functions v{t) and w{t) absolutely continuous on 
[0,T], one has the equality: 



t 5 « 

f3 f f v'{r])dr] f w'{s)ds 



r{i-/3)J {t-O'-^J {t-r,YJ {t-sY' 





(4) 



where < /3 < 1. 

Proof. Let us consider the difference 



1 f w'{s){v{t) - v{s)) + v'{s){w{t) - w{s)) 



it-sY 



ds 



t t 
1 f v'{s)ds 



T{l-f3)J {t-sY 





t t 

1 /" w'{s)ds 



r(i-/3)i {t-sY 





1 /" ,/ .X /" v'{s)ds 



+ 



t « 
1 /" ,/ .X /" w'{s)ds 



r(i-/3)7 V {t-sY r(i-/3)7 V {t-sy 
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v'{rf)d'q f w'{s)ds 



1 f(t-eY^ / 

\o , 

t ? « 

(3 f d^ f v'{r])dr] f w'{s)ds 



d^ 



r(i-/3)y (t-o^-^y (t-r^)^; (t-s)^' 



The proof of the lemma is complete. 

If v{t) = w{t) then from the Lemma 1 one has the following: 
Corollary 1. For any function v{t) absolutely continuous on [0,T], the 

following equality takes place: 

\o / 

where < (3 < 1. 

I 

Let us use the following notation: ||m||q = J u'^{x,t)dx, DQfu{x,t) = 



t 

j {t — s)^^^-u(x, s)ds/r{l3) - fractional Riemann-Liouville integral of order /3. 



Theorem 1. li k{x,t) E C^'%Qt), q{x,t), f{x,t) e C{Qt), k{x,t) > 
Ci > 0, q{x, t) >0 everywhere on Qt, then the solution u{x, t) of the problem 
©^dS]) satisfies the a priori estimate: 



Dq/^^ ^u^{x,t)dx + Ci / s) llgfis < 





t I 



Proof. Let us multiply equation ([1]) by u{x, t) and integrate the resulting 
relation over x from to /: 
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+ J q{x,t)u {x,t)dx = J u{x,t)f{x,t)dx. 



Then transform the terms in identity ([7]) as: 



(7) 







u{x, t)f{x, t)dx 



<e\\u{xMl + ^J\f{^Ml ^>0, 



Using the equahty ([5]) one obtains 



u{x,t)dQl;^\{x,t)dx > - I Bq^'^' 'u^{x,t)dx. 



Taking into account the above performed transformations, from the identity 
(I7j) one arrives at the inequahty 



\ 1 9otV(a:,t)da: + Ci||«.(x,t)||^<5||n(a:,t)||^ + l||/(x,t)||g. (8) 



Using the inequahty ||M(x,t)||Q < (/^/2)||u2,.(x, t)||Q, from the inequahty 
at e = Ci/l^ one obtains 



J d^t^u\x,t)dx + cr\\u.{x,ml<^y{^Ml 



(9) 



Changing the variable t by s in the inequahty (j9]) and integrating it over 
s from to t, one obtains a priori estimate ([6]). The uniqueness and the 
continuous dependence of the solution of the problem ([TJ-Q on the input 
data follow from the a priori estimate (EI). 



5 



The solution of the problem ([I])-© with a{x) = a {a = const) satisfies 
the a priori estimates: 



t t 
D'^-^\\u{x,t)\\l+ci J \Mx,s)\\lds < — J \\fix,s)\\lds+Y^j-—^\\uoix)\\l, 





(10) 



\Hx,t)\\l + D,,^u,{x,t)\\l < M {D^,^fix,t)\\l + \\ul{x)\\l) . (11) 

Inequality ( ITOl) follows from ([6]), and the a priori estimate (fTTi) follows 
from inequality (E]) with a{x) = a. Actually, applying the fractional inte- 
gration operator Dq^°' to the both sides of inequality ([8]), one arrives at the 
estimate (ITTl) . which contains the constant M = max{/^/ci, 1}/ min{l,ci}. 

2.2. The Robin boundary value problem. 

In the problem ([I])-(l3]) we replace the boundary conditions with 

r k{0,t)u,{0,t) = P,{t)u{0,t)-fi,{t), 

\ ~k{l, t)u^{l, t) = P2{tHl, t) - /X2(t). 

In the rectangle Qt we consider the Robin boundary value problem ([1]), 

(ED, m- 

Theorem 2. If k{x,t) G C^'°{Qt), q{x,t), f{x,t) G C{Qt), k{x,t) > 
ci > 0, q{x,t) > everywhere on Qt, (3i{t),Hi{t) G C[0,T], (3i{t) > /3o > 0, 
for all t G [0,T], i = 1,2, then the solution u{x,t) of the problem ([1]), ([3]), 
(fl2|) satisfies the a priori estimate: 



D°/")-V(x,t)da; + 7 / (||M,.(x,s)||2 + n2(0,s) + n2(/,s))rfs 



< 



5 

< - 

7 



(||/(x,s)||2 + /i2(s)+/i2(s))rfs 



V{2-a{x)) 



u1{x)dx, (13) 
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where 7 = min{ci, /3o}, S = max{l + /, l^}. 

Proof. Multiply the equation ([T]) by u{x,t) and integrate the resulting 
relation over x from to /: 



(14) 



DO 

Then, transform the terms of the identity f|T4l) : 

J u{x,t)dQ}^^u{x,t)dx > - J dQ^'^^u'^{x,t)dx. 



- / u{ku:,)^dx = i3i{t)u'^{0,t)+(32{t)u'^{l,t)-iJ,i{t)u{0,t)-iJ,2it)u{l,t)+ / kuldx, 



ufdx 



- ^\m\o + -^wj lie 



e>0. 



From f|T4|) . taking into account the transformations performed, one arrives 
at the inequality 



9o/"^m'(x, t)dx + c4u,{x, t) \\l + Pou'{0, t) + Pou\l, t) < 



< /ii(t)n(0, t) + fi2it)uil, t) + e\\u\\1 + ^WfWl (15) 

Using the inequalities fii{t)u{0,t) < 6u'^{0,t) + (4e)"Vi(^)? f^2{t)u{l,t) < 
eu^{l,t) + (4e)-V2W> ^ > 0; lk(a;,i)|lo < l'^\\u^{x,i)\\l + l{u'^{Q,t) + u^il.t)) 
with e = 7/(25), from (fT5|) one has the following inequality 

I 

do^^'^u^ix, t)dx + 7 (||mx(x, t) 11^ + ^2(0, t) + u\l, t)) < 







<^(||/(x,t)||^ + /.?(t)+/.^(t)). (16) 



Changing variable t by s in inequality f ll6p and integrating it over s from 
to t, one obtains the a priori estimate ( IT^ . 

The uniqueness and the continuous dependence of the solution of problem 
dl]), ©, f|T2l) on the input data follow from the a priori estimate f|T3|) . 

The solution of the problem ([1]), ([3]), f|T2l) with o;(a;) = a (a = const) 
satisfies the following a priori estimates: 



t 




D^i'\\u{x,t)\\l + l\^j {\\u,{x,s)\\l + u\0,s)+u\l,s))dsj < 

(ll/(^,^)llo + /^?(^)+/^2(^))rfn +r|^lko(x)|^^ (17) 

/ 

t)\\l + DoT llo < t) 
+ D,,yiit) + D,,ylit) + \\ulix)\\l), (18) 

Inequality (fT7|) follows from (1131) . and the a priori estimate (1181) follows from 
inequality (ITB]) with a(x) = a. Actually, applying the fractional integration 
operator D^^"' to the both sides of inequality (ITB]) . one arrives at the estimate 
f[T8|) . which contains the constant M = max{5/7, 1}/ min{l,7}. 

3. Boundary value problems in difference setting 

3.1. The Dirichlet boundary value problem 

Suppose that a solution u{x,t) G C'^'^{Qt) of the problem ([I])-© exists, 
and the coefficients of the equation ([1]) and the functions f{x, t), uo{x) satisfy 
the smoothness conditions, required for the construction of difference schemes 
with the order of approximation 0{t + h"^). 

In the rectangle Qt we introduce the grid Uhr = oJh ^ ^t, where Uh = 

= ih,i = 0, 1, A^, hN = I}, Ur = {tj =jT,j = 0, 1, jo,rjo = T}. 

To problem ([I])-(l3]), we assign the difference scheme: 

A^i^y = A{ay^+' + {l-a)y^)+^, i = l,...,N-l, j = 1, jo - 1, (19) 
2/(0,t) = 0, 2/(/,t)=0, j = 0,...,jo, (20) 
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y{x,0) = uo{x), i = 0,...,N, (21) 

where Ay = {ays)x-dy, a = A;(xi_i/2, t), d = q{xi,i), ip = f{xi,i), i = tj+1/2, 
< a < 1, A^l^y = Ei=oit'P+i - t]ZT)yt/^i'^ - a.) - a diference analogue 
of the Caputo fractional derivative of order ) UM. 

According to [l^, [2^ the difference scheme fll9l) - fl21l) has the order of 
approximation 0{t + h'^). 

Lemma 2. For every function y{t) defined on the grid COr one has the 
inequahties 

y^+'A-,,y > (y^) + ^^^^^|^(A-|/)^ (22) 

y^A-,,y > ^-A^y') - ^^^i^T-") (^otl/)^- (23) 
Proof. Inequahty (!22|) is equivalent to the inequality 



r(2-Q)^ 2 2 



s — — 1_ s — 

- i!Q|^(a;„)^ > 0. (24) 

Here we consider the sums to be equal to zero if the upper summation 
limit is less than the lower one. 
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Let us introduce the following notation: Yl'l=o(^]-s+i ~ '^]~s)yt — 'V^'^^, 
tlienyO = (tj;?-tj-")-V, = Tit]zt+, - t]zty'vt A: = 1, 2, j. 
Taking into account the introduced notation, we rewrite the inequality ([24 
as 

j 



r(2-a)^ ^^-"^^ * 2r(2-a; 



21(2 -«)^ ^ 21(2- a) 

-(i]:."-i]Xi)"^)(^'^T>0. (25) 

Obviously, inequality ([25]) is valid since (t]:fc+i-t]:fc)"^-(t]Zfc -tjl^i)"^ > 
0,fc = 0,l,...,j-l. 

Let us prove now the inequality fl25]) . Since = y^^^ — ryt, one obtains 

y^A-,,y - ^-A-Uy') + ^^^^{A^oty? = 
= y^^'Aly - -Al{y') + ^J^ {A'^,,yf - ry.A^y = 

2r(2-a)(2-2i-)^ ' r(2-«) * ^2r(2 - «) ^ ^-^=+1 ^-^^ 
y^j-k ^j-k-i) ){v ) - 9rr9_^V9_9i-"^^^ > ^r(o^^\^ ^ + 



2r(2-a)(2-2i-°)' ' r(2-a) 

T°(2-2^ ") + I Vrfft^-" -t^-^r^- 

2r(2-a)(2i-°-l)^^ ^2r(2-«)^^^^-'=+i ^-'^^ 



fc=0 
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The proof of the lemma is complete. 

Corollary 2. For any function y{t) defined on the grid a)^ one has the 
inequality 

1 T°'V(2 — rv) 

{cry^^' + (1 - <j)y' W,ty > 2^0^(1/') + _ 21-^) - 2^"")^ " l) i^otV)'- 

(26) 

Theorem 3. The difference scheme ([T9])-(l21]) at a > 1/(3 - 2i-'°'^("«)) 
is absolutely stable and its solution satisfies the following a priori estimate: 

(f(2^'B^-?--^]:;^0(/-)^^ 

+Ci^||ayi'+^ + (l-a)i/|']|^r< 

j'=0 

where {y,v) = Ya~i^ ViVih, {y,v\ = ^^^yiVih, \\y\\l = {y,y), \\y]\l = 

{y,y]- 

Proof. Let us multiply scalarly equation ( |T9l) by y*-'^^ = ay^^^ + (l — a)y^: 
(AS,.2/,2/('^))-(Ay('^\y('^)) = (^,l/('^)). (28) 
Let us transform the terms in identity (!28l) : 



(l/^'^^AS,l/)>^(l,AS,.(^/'))+ 

+( ;(2^!:j ((3-2^-0^-l),(Ag^y)^ 
-(Ay('^),y(-)) = (a,(yf )2] + (ci, (y('^))2), 



Taking into account the above-performed transformations, from identity (l28ll 
at (T > 1/(3 — 2-'^^™™^"')) one arrives at the inequality 

^(1, AS^.(y^)) + cilbf ]|^ < eWy^^^^Wl + j-JMl (29) 
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From (^Uij at e = Ci//^, using that \\y\\l < {l"^ /2)\\ys]\l, one obtains the 
inequahty 

(l,A-(y2))+Ci|bf]|^<|-||v.||^. (30) 

Multiplying the inequality ( l30l) by r and summing over j' from to j, 
one obtains the a priori estimate (|2711 . The stability and convergence of the 
difference scheme ( !T9|) - fl2Tl) follow from the a priori estimate fl27|) . 

If a{x) = a {a = const) then the solution of the problem f|T9l) - fl2Tl) 
satisfies the following a priori estimate: 



i'=o i'=o 

i j'=0 ^ ^ 

Here the results are obtained for the homogeneous boundary conditions 
u{0,t) = 0, u{l,t) = 0. In the case of inhomogeneous boundary conditions 
u{0,t) = /ii(t), u{l,t) = fi2it) the boundary conditions for the difference 
problem will have the following form: 

yiO,t) = fiiit), yil,t) = fi2{t). (32) 



Convergence of the difference scheme ( ]T9l) . (I2T1) . (1321) follows from the 
results obtained above. Actually, let us introduce the notation y = z + u. 
Then the error z = y — u is a solution of the following problem: 

A-Sl^z = A{az^+' + {l-cx)z^)+^, t = l,...,N-l, j = 1, jo - 1, (33) 

z(0,t) = 0, zil,t) = 0, j = 0,...,jo, (34) 

z{x,0) = 0, i = 0,...,N, (35) 

where ip = A{au^+'^ + (1 - (^)u^) - ^ot + V = 0(r + h"^). 

The solution of the problem f p3|) -f p5|) satisfies the estimation f l27|) so that 
the solution of the difference scheme f lT9|) . f l2T]) . fl32|) converges to the solution 
of the corresponding differential problem with order 0(r + /i^). 
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3.2. Numerical results 

In this section, the following variable order time fractional diffusion equa- 
tion is considered: 



where Q!(a;)— -^^ , k{x,f) — (^^,,£2^iu^,,^A^2x+i]+(x^+x+i)(i2x-i+2)^ — 



Uo{x) = {x^ + X + l)(3a;^ + 2x + 1). 

The exact solution is u(x, t) ^ (x^ + x + l){3x^ + 2x + l){t^ + + 1). 

All the calculations are performed at cr = 1/(3 — 2-'^"™'"'^"')), where for the 
considered example min(Q;j) = 1/10. 

A comparison of the numerical solution and exact solution is provided in 



Table 2 shows that if h — 1/500, then as the number of time of our 
approximate scheme is decreased, a redaction in the maximum error takes 
place, as expected and the convergence order of time is 0(r), where the 
convergence order is given by the following formula: Convergence order = 

Table 3 shows that when we take h'^ = r, as the number as spatial subin- 
tervals/time steps is decreased, a reduction in the maximum error takes place, 
as expected the convergence order of the approximate scheme is 0{h^ + r), 
where the convergence order is given by the following formula: Convergence 
order = log^i — . 




(36) 




Table 1. 
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Table 1 



The error, numerical solution and exact solution, when t = 1, h 


= 1/10, r = 1/100. 


bpace(Xij 


Numerical solution 


Exact solution 


Error 


0.0000 


5.0000000 


5.0000000 


0.0000000 


0.1000 


6.6068921 


6.6076515 


0.0007594 


0.2000 


8.4813103 


8.4849920 


0.0036817 


U.oUUU 


iU. / /ooy 


1 n 7770Qnc^ 

iU. i i 1 ZoUO 


0.0094736 


n A nnn 
U.4UUU 


1 Q 71 m o/i 
ic). ( iy iy4o 


1 Q 7Qfi1 7Rn 


n ni SOS1 K 
U.UioyoiO 


n c;nnn 
U.JdUUU 


LI . 14:0 lOOO 


1 7 77Q /I Q7C^ 


u.uozoouy 


O.dOOO 


23.4561082 


oo rn/^on/in 

23.5063040 


U.UOUiyOo 


n 7nnn 
U. ( uuu 


Oi.oU4i ( /D 


oi .o ( 00040 


n ncn/:?nin 
0.0696919 


0.8000 


44.1761384 


44.2609280 


u.uo4 ( oyu 


n nnnn 

u.yuuu 


CO <^ni fi7c;Q 


o/.o / yoUoo 


0.0776277 


1 nnnn 


90.0000000 


90.0000000 


n nnnnnnn 
u.uuuuuuu 










Maximum error 


behavior versus time grid size 


reduction at t = 


1 when h = 1/500. 


T 


Maximum error 




Convergence order 


1/256 


0.0344960 






1/1024 


0.0086690 




0.996 


1/4096 


0.0021738 




0.998 


Table 3 








Maximum error 


behavior versus grid size reduction at t = 1 when = r. 


h 


Maximum error 




Convergence order 


1/40 


0.0056275 






1/80 


0.0014141 




1.993 


1/160 


0.0003542 




1.997 



3.3. The Robin boundary value problem. 

To the differential problem ([I]) , (EJ , (fT2|) we assign the following difference 
scheme: 

A^y = A{ay^+' + (1 - a)yJ) + ^, t = 0, N,j = 1, jq, (37) 
y{x,0) = Uo{x), t = 0,...,N, (38) 
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where Ky = {aiy^-Piy)/ {0.5h),i = 0, Ay = {ays)x-dy, i = 1, N-1, 
Ay = {ttNyx - §2y)/{0.5h),i = N, Lpo = {2fli)/h, LfN = {2fl2)/h, /3i = 
/9i + O.bhdo, 132 = (32 + 0.5hdN, Jxi = + 0.5/i/o, /i2 = + 0.5hfN- 
The difference scheme fl371) - fl38l) has the order of approximation 0{t + h'^). 

Theorem 4. The difference scheme ([37])-([38l) at a > 1/(3 - 2^-"'"'^''^^) 
is absolutely stable and its solution satisfies the following a priori estimate: 



1 



+ 



j'=0 



M^j'+1^2 



< 



T+ 



^ i'=0 



+ 



r(2 - ai 



(39) 



where 7 = min{ci, /3o}, ^ = max{l+/, f}, [y, v] = Y^iLi^ yiVih+0.5yoVoh+ 
O.bvNVNh, \[y]\l = [y,y], {y'^'^Y^' = ^V''^' + (1 " ^)/- 

Proof. Let us multiply scalarly equation (1371) by y^'^^ = cry^^^ + {l — a)y^: 



Let us transform the terms occurring in identity ( HOj) 



(40) 



+ ((3-2^""V-l).(Ag^yn, 
-[Al/(-),y(-)] = (3^y', + hvl + (a, (z/i"^^)'] + [rf, d/^'^^)'], 

\\^>.y^''\ < eU'^X + /ii2/o + l^2yN + ^ll^ll^, e > 0. 

Taking into account the above performed transformations, from identity (140!) 
at (7 > 1/(3 — 2-'^^™™("»)) one arrives at the inequality 



2[l'A^(^')] + cilbf]|^ + /3o(yo' + y^)< 
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< ^(Ib^'^^llo + vo + yl) + ^(-^? + Ml)- (41) 

From fj40l) aX e = 7/(25), using that \\y\\l < /^||ys]lo + Uz/o + ^"^^ 
the following inequality: 

[l,A-(y2)]+7(||l/f]|^ + l/o' + 2/^)< 

<-(/i?+/i2+ll</'llo)- (42) 
7 

Multiplying inequality fH2|) by r and summing over j' from to j, one obtains 
a priori estimate (!39l) . The stability and convergence of the difference scheme 
(!371) -(l38l) follow from the a priori estimate ( |39ll . 

If a(a;) = a (a = const), then for the solution of the problem (|371) - (!38ll 
one has the following a priori estimate: 

+7 E {\\i^y'^'y'^']\l + iiy'oY^'r + ((?/{.^^)^"^^)0 - < 

i'=o 

<-E((/^^^'')' + (/^f'''')' + ii^^'iio)^ 



^ i'=o 



+ ^;|l^lK(x.)]|g. (43) 



4. Conclusion 

The results obtained in the present paper allow to apply the method of 
energy inequalities to finding a priori estimates for boundary value problems 
for the fractional diffusion equation in differential and difference settings 
exactly as in the classical case {a{x) = 1). It is interesting to note that the 
condition a > 1/(3 - 2^""') at a{x) = 1 turns into the well known condition 
0" > 1/2 of the absolute stability of the difference schemes for the classical 
diffusion equation. 
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